arXiv:1506.07834v2 [nucl-ex] 24 Jan 2016 


Measurement of higher cumulants of net-charge multiplicity distributions in Au+Au 

collisions at y' s NN = 7.7—200 GeV 

A. Adare, 13 S. Afanasiev, 31 C. Aidala, 40,44,45 N.N. Ajitanand, 65 Y. Akiba, 59,60 R. Akimoto, 12 H. Al-Bataineh, 53 

J. Alexander, 65 H. Al-Ta’ani, 53 A. Angerami, 14 K. Aoki, 33,36,59 N. Apadula, 29,66 Y. Aramaki, 12,59 H. Asano, 36,59 

E. C. Aschenauer, 7 E.T. Atomssa, 37,66 R. Averbeck, 66 T.C. Awes, 55 B. Azmoun, 7 V. Babintsev, 25 M. Bai, 6 

G. Baksay, 20 L. Baksay, 20 B. Bannier, 66 K.N. Barish, 8 B. Bassalleck, 52 A.T. Basye, 1 S. Bathe, 5 ’ 8 ’ 60 
V. Baublis, 58 C. Baumann, 46 S. Baumgart, 59 A. Bazilevsky, 7 S. Belikov, 7, 0 R. Belmont, 13 ’ 45 ’ 70 R. Bennett, 66 

A. Berdnikov, 62 Y. Berdnikov, 62 A.A. Bickley, 13 D. Black, 8 D.S. Blau, 35 J.S. Bok, 52 ’ 53 ’ 74 K. Boyle, 60 ’ 66 
M.L. Brooks, 40 J. Bryslawskyj, 5 H. Buesching, 7 V. Bumazhnov, 25 G. Bunce, 7,60 S. Butsyk, 40,52 C.M. Camacho, 40 

S. Campbell, 14 ’ 66 P. Castera, 66 C.-H. Chen, 60 ’ 66 C.Y. Chi, 14 M. Chiu, 7 I.J. Choi, 26 ’ 74 J.B. Choi, 10 S. Choi, 64 

R.K. Choudhury 4 P. Christiansen, 42 T. Chujo, 69 P. Chung, 65 O. Chvala, 8 V. Cianciolo, 55 Z. Citron, 66,72 

B. A. Cole, 14 M. Connors, 66 P. Constantin, 40 N. Cronin, 47,66 N. Crossette, 47 M. Csanad, 18 T. Csorgo, 73 
T. Dahms, 66 S. Dairaku, 36,59 I. Danchev, 70 K. Das, 21 A. Datta, 44,52 M.S. Daugherity, 1 G. David, 7 

K. Dehmelt, 20,66 A. Denisov, 25 A. Deshpande, 60,66 E.J. Desmond, 7 K.V. Dharmawardane, 53 O. Dietzsch, 63 

L. Ding, 29 A. Dion, 29,66 J.H. Do, 74 M. Donadelli, 63 L. D’Orazio, 43 O. Drapier, 37 A. Drees, 66 K.A. Drees, 6 

J. M. Durham, 40,66 A. Durum, 25 D. Dutta, 4 S. Edwards, 6,21 Y.V. Efremenko, 55 F. Ellinghaus, 13 T. Engelmore, 14 

A. Enokizono, 39 ’ 55,59,61 H. En’yo, 59,60 S. Esurni, 69 K.O. Evser, 7,8 B. Fadem, 47 D.E. Fields, 52 M. Finger, 9 

M. Finger, Jr., 9 F. Fleuret, 37 S.L. Fokin, 35 Z. Fraenkel, 72 ’|j J.E. Frantz, 54,66 A. Franz,' A.D. Frawley, 21 

K. Fujiwara, 59 Y. Fukao, 59 T. Fusayasu, 49 K. Gainey, 1 C. Gal, 66 P. Garg, 3 A. Garishvili, 67 I. Garishvili, 39,67 

F. Giordano, 26 A. Glenn, 13,39 H. Gong, 66 X. Gong, 65 M. Gonin, 37 Y. Goto, 59,60 R. Granier de Cassagnac, 37 

N. Grau, 2,14 S.V. Greene, 70 M. Grosse Perdekamp, 26,60 Y. Gu, 65 T. Gunji, 12 L. Guo 40 H.-A. Gustafsson , 42, 0 

T. Hachiya, 24,59 J.S. Haggerty, 7 K.I. Hahn, 19 H. Hamagaki, 12 J. Hamblen, 67 R. Han, 57 J. Hanks, 14,66 
E.P. Hartouni, 39 K. Hashimoto, 59,61 E. Haslum, 42 R. Hayano, 12 S. Hayashi, 12 X. He, 22 M. Heffner, 39 
T.K. Hemmick, 66 T. Hester, 8 J.C. Hill, 29 M. Hohlmann, 20 R.S. Hollis, 8 W. Holzmann, 14 K. Homma, 24 B. Hong, 34 
T. Horaguchi, 24,69 Y. Hori, 12 D. Hornback, 67 S. Huang, 70 T. Ichihara, 59,60 R, Ichimiya, 59 J. Ide, 47 H. Iinuma, 33 
Y. Ikeda, 59,69 K. Irnai, 30,36,59 Y. Imazu, 59 J. Imrek, 17 M. Inaba, 69 A. Iordanova, 8 D. Isenhower, 1 M. Ishihara, 59 
A. Isinhue, 47 T. Isobe, 12,59 M. Issah, 70 A. Isupov, 31 D. Ivanishchev, 58 B.V. Jacak, 66 M. Javani, 22 J. Jia, 7,65 
X. Jiang, 40 J. Jin, 14 B.M. Johnson, 7 K.S. Joo, 48 D. Jouan, 56 D.S. Jumper, 1,26 F. Kajihara, 12 S. Kametani, 59 
N. Kamihara, 60 J. Kamin, 66 S. Kaneti, 66 B.H. Kang, 23 J.H. Kang, 74 J.S. Kang, 23 J. Kapustinsky 40 

K. Karatsu, 36,59 M. Kasai, 59,61 D. Kawall 44,60 M. Kawashima, 59,61 A.V. Kazantsev, 35 T. Kempel, 29 J.A. Key, 52 
P.K. Khandai, 3 A. Khanzadeev, 58 K.M. Kijima, 24 B.I. Kim, 34 C. Kim, 34 D.H. Kim, 48 D.J. Kim, 32 E. Kim, 64 

E.-J. Kim, 10 H.J. Kim, 74 K.-B. Kim, 10 S.H. Kim, 74 Y.-J. Kim, 26 Y.K. Kim, 23 E. Kinney, 13 K. Kiriluk, 13 
A. Kiss, 18 E. Kistenev, 7 J. Klatsky, 21 D. Kleinjan, 8 P. Kline, 66 L. Kochenda, 58 Y. Komatsu, 12,33 B. Kornkov, 58 
M. Konno, 69 J. Koster, 26,60 D. Kotchetkov, 52,54 D. Kotov, 58,62 A. Kozlov, 72 A. Krai, 15 A. Kravitz, 14 F. Krizek, 32 

G.J. Kunde, 40 K. Kurita, 59,61 M. Kurosawa, 59,60 Y. Kwon, 74 G.S. Kyle, 53 R. Lacey, 65 Y.S. Lai, 14 J.G. Lajoie, 29 
A. Lebedev, 29 B. Lee, 23 D.M. Lee, 40 J. Lee, 19 K. Lee, 64 K.B. Lee, 34,40 K.S. Lee, 34 S.H. Lee, 66 S.R. Lee, 10 
M.J. Leitch 40 M.A.L. Leite, 63 M. Leitgab, 26 E. Leitner, 70 B. Lenzi, 63 B. Lewis, 66 X. Li, 11 P. Liebing, 60 S.H. Lim, 74 

L.A. Linden Levy, 13,39 T. Liska, 15 A. Litvinenko, 31 H. Liu, 40,53 M.X. Liu, 40 B. Love, 70 R. Luechtenborg 46 
D. Lynch, 7 C.F. Maguire, 70 Y.I. Makdisi, 6 M. Makek, 72,75 A. Malakhov, 31 M.D. Malik, 52 A. Manion, 66 
V.I. Manko, 35 E. Mannel, 7,14 Y. Mao, 57,59 T. Maruyama, 30 H. Masui, 69 S. Masumoto, 12,33 F. Matathias, 14 

M. McCumber, 13 ’ 40,66 P.L. McGaughey, 40 D. McGlinchey, 13,21 C. McKinney, 26 N. Means, 66 A. Meles, 53 
M. Mendoza, 8 B. Meredith, 26 Y. Miake, 69 T. Mibe, 33 J. Midori, 24 A.C. Mignerey 43 P. Mikes, 9,28 K. Miki, 59,69 

A. Milov, 7,72 D.K. Mishra, 4 M. Mishra, 3 J.T. Mitchell, 7 Y. Miyachi, 59,68 S. Miyasaka, 59,68 A.K. Mohanty, 4 

S. Mohapatra, 65 H.J. Moon, 48 Y. Morino, 12 A. Morreale, 8 D.P. Morrison,' '0 M. Moskowitz, 47 S. Motschwiller, 47 

T.V. Moukhanova, 35 T. Murakami, 36,59 J. Murata, 59,61 A. Mwai, 65 T. Nagae, 36 S. Nagamiya, 33,59 J.L. Nagle, 13, @ 
M. Naglis, 72 M.L Nagy, 18,73 I. Nakagawa, 59,60 Y. Nakamiya, 24 K.R. Nakamura, 36,59 T. Nakamura, 33,59 
K. Nakano, 59,68 C. Nattrass, 67 A. Nederlof, 47 P.K. Netrakanti, 4 J. Newby, 39 M. Nguyen, 66 M. Nihashi, 24,59 
T. Niida, 69 R. Nouicer, 7,60 N. Novitzky, 32,66 A. Nukariya, 12 A.S. Nyanin, 35 H. Obayashi, 24 E. O’Brien, 7 
S.X. Oda, 12 C.A. Ogilvie, 29 M. Oka, 69 K. Okada, 60 Y. Onuki, 59 A. Oskarsson, 42 M. Ouchida, 24,59 K. Ozawa, 12,33 
R. Pak, 7 V. Pantuev, 27,66 V. Papavassiliou, 53 B.H. Park, 23 I.H. Park, 19 J. Park, 10,64 S. Park, 64 S.K. Park, 34 



2 


W.J. Park, 34 S.F. Pate, 53 L. Patel, 22 H. Pei, 29 J.-C. Peng, 26 H. Pereira, 16 D.V. Perepelitsa, 7,14 V. Peresedov, 31 
D.Yu. Peressounko, 35 R. Petti, 7,66 C. Pinkenburg, 7 R.P. Pisani, 7 M. Proissl, 66 M.L. Purschke, 7 A.K. Purwar, 40 
H. Qu, 1,22 J. Rak, 32 A. Rakotozafindrabe, 37 I. Ravinovich, 72 K.F. Read, 55,67 K. Reygers 46 D. Reynolds, 65 
V. Riabov, 51,58 Y. Riabov, 58,62 E. Richardson, 43 N. Riveli, 54 D. Roach, 70 G. Roche, 41 B S.D. Rolnick, 8 
M. Rosati, 29 C.A. Rosen, 13 S.S.E. Rosendahl, 42 P. Rosnet, 41 P. Rukoyatkin, 31 P. Ruzicka, 28 M.S. Ryu, 23 
B. Sahlmueller, 46,66 N. Saito, 33 T. Sakaguchi, 7 K. Sakashita, 59,68 H. Sako, 30 V. Samsonov, 51,58 M. Sano, 69 
S. Sano, 12,71 M. Sarsour, 22 S. Sato, 30,33 T. Sato, 69 S. Sawada, 33 K. Sedgwick, 8 J. Seele, 13 R. Seidl, 26,59,60 
A.Yu. Semenov, 29 A. Sen, 22,67 R. Seto, 8 P. Sett, 4 D. Sharma, 66,72 I. Shein, 25 T.-A. Shibata, 59,68 K. Shigaki, 24 
M. Shimomura, 29,50,69 K. Slioji, 36,59 P. Shukla, 4 A. Sickles, 7,26 C.L. Silva, 29,40,63 D. Silvermyr, 42,55 C. Silvestre, 16 
K.S. Sim, 34 B.K. Singh, 3 C.P. Singh, 3 V. Singh, 3 M. Skolnik, 47 M. Slunecka, 9 S. Solano, 47 R.A. Soltz, 39 
W.E. Sondheim, 40 S.P. Sorensen, 67 I.V. Sourikova, 7 N.A. Sparks, 1 P.W. Stankus, 55 P. Steinberg, 7 E. Stenlund, 42 
M. Stepanov 44,53, B A. Ster, 73 S.P. Stoll, 7 T. Sugitate, 24 A. Sukhanov, 7 J. Sun, 66 J. Sziklai, 73 E.M. Takagui, 63 
A. Takahara, 12 A. Taketani, 59,60 R. Tanabe, 69 Y. Tanaka, 49 S. Taneja, 66 K. Tanida, 36,59,60,64 M.J. Tannenbaum, 7 
S. Tarafdar, 3,72 A. Taranenko, 51,65 P. Tarjan, 17 E. Tennant, 53 H. Tliemann, 66 T.L. Thomas, 52 T. Todoroki, 59,69 
M. Togawa, 36,59 A. Toia, 66 L. Tomasek, 28 M. Tomasek, 15,28 H. Torii, 24 R.S. Towell, 1 I. Tserruya, 72 
Y. Tsuchimoto, 12,24 T. Tsuji, 12 C. Vale, 7,29 H. Valle, 70 H.W. van Hecke, 40 M. Vargyas, 18 E. Vazquez-Zambrano, 14 
A. Veicht, 14,26 J. Velkovska, 70 R. Vertesi, 17,73 A.A. Vinogradov, 35 M. Virius, 15 B. Voas, 29 A. Vossen, 26 
V. Vrba, 15,28 E. Vznuzdaev, 58 X.R. Wang, 53,60 D. Watanabe, 24 K. Watanabe, 59,61,69 Y. Watanabe, 59,60 
Y.S. Watanabe, 12 F. Wei, 29,53 R. Wei, 65 J. Wessels, 46 S. Whitaker, 29 S.N. White, 7 D. Winter, 14 S. Wolin, 26 
J.P. Wood, 1 C.L. Woody, 7 R.M. Wright, 1 M. Wysocki, 13,55 B. Xia, 54 W. Xie, 60 Y.L. Yamaguchi, 12,59,66 
K. Yamaura, 24 R. Yang, 26 A. Yanovich, 25 J. Ying, 22 S. Yokkaichi, 59,60 Z. You, 40,57 G.R. Young, 55 
I. Younus, 38,52 I.E. Yushmanov, 35 W.A. Zajc, 14 A. Zelenski, 6 C. Zhang, 55 S. Zhou, 11 and L. Zolin 31 

(PHENIX Collaboration) 

1 Abilene Christian University, Abilene, Texas 79699, USA 
2 Department of Physics, Augustana University, Sioux Falls, South Dakota 57197, USA 
3 Department of Physics, Banaras Hindu University, Varanasi 221005, India 
A Bhabha Atomic Research Centre, Bombay fOO 085, India 
5 Baruch College, City University of New York, New York, New York, 10010 USA 
6 Collider-Accelerator Department, Brookhaven National Laboratory, Upton, New York 11973-5000, USA 
'Physics Department, Brookhaven National Laboratory, Upton, New York 11973-5000, USA 
8 University of Califomia-Riverside, Riverside, California 92521, USA 
9 Charles University, Ovocny trh 5, Praha 1, 116 36, Prague, Czech Republic 
10 Chonbuk National University, Jeonju, 561-756, Korea 
11 Science and Technology on Nuclear Data Laboratory, China Institute 
of Atomic Energy, Beijing 102413, People’s Republic of China 
12 Center for Nuclear Study, Graduate School of Science, University of Tokyo, 7-3-1 Hongo, Bunkyo, Tokyo 113-0033, Japan 

13 University of Colorado, Boulder, Colorado 80309, USA 
14 Columbia University, New York, New York 10027 and Nevis Laboratories, Irvington, New York 10533, USA 
15 Czech Technical University, Zikova j, 166 36 Prague 6, Czech Republic 
16 Dapnia, CEA Saclay, F-91191, Gif-sur-Yvette, France 
17 Debrecen University, H-jOlO Debrecen, Egyetem ter 1, Hungary 
1S ELTE, Eotvos Lorand University, H-1117 Budapest, Pazmany P. s. 1/A, Hungary 
19 Ewha Womans University, Seoul 120-750, Korea 
20 Florida Institute of Technology, Melbourne, Florida 32901, USA 
21 Florida State University, Tallahassee, Florida 32306, USA 
22 Georgia State University, Atlanta, Georgia 30303, USA 
23 Hanyang University, Seoul 133-792, Korea 
24 Hiroshima University, Kagamiyama, Higashi-Hiroshima 739-8526, Japan 
25 IHEP Protvino, State Research Center of Russian Federation, Institute for High Energy Physics, Protvino, If2281, Russia 
26 University of Illinois at Urbana- Champaign, Urbana, Illinois 61801, USA 
27 Institute for Nuclear Research of the Russian Academy of Sciences, prospekt 60-letiya Oktyabrya 7a, Moscow 117312, Russia 
28 Institute of Physics, Academy of Sciences of the Czech Republic, Na Slovance 2, 182 21 Prague 8, Czech Republic 

29 Iowa State University, Ames, Iowa 50011, USA 
30 Advanced Science Research Center, Japan Atomic Energy Agency, 2-4 
Shirakata Shirane, Tokai-mura, Naka-gun, Ibaraki-ken 319-1195, Japan 
31 Joint Institute for Nuclear Research, 141980 Dubna, Moscow Region, Russia 
32 Helsinki Institute of Physics and University of Jyvaskyla, P.O.Box 35, FI-40014 Jyvaskyla, Finland 
33 KEK, High Energy Accelerator Research Organization, Tsukuba, Ibaraki 305-0801, Japan 


3 


34 Korea University, Seoul, 136-701, Korea 
3o National Research Center “Kurchatov Institute”, Moscow, 123098 Russia 
36 Kyoto University, Kyoto 606-8502, Japan 

3l Laboratoire Leprince-Ringuet, Ecole Polytechnique, CNRS-IN2P3, Route de Saclay, F-91128, Palaiseau, France 
38 Physics Department, Lahore University of Management Sciences, Lahore 54792, Pakistan 
39 Lawrence Livermore National Laboratory, Livermore, California 94550, USA 
40 Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA 
41 LPC, Universite Blaise Pascal, CNRS-IN2P3, Clermont-Fd, 63177 Aubiere Cedex, France 
42 Department of Physics, Lund University, Box 118, SE-221 00 Lund, Sweden 
43 University of Maryland, College Park, Maryland 20742, USA 
44 Department of Physics, University of Massachusetts, Amherst, Massachusetts 01003-9337, USA 
45 Department of Physics, University of Michigan, Ann Arbor, Michigan 48109-1040, USA 
46 Institut fur Kernphysik, University of Muenster, D-48149 Muenster, Germany 
47 Muhlenberg College, Allentown, Pennsylvania 18104-5586, USA 
48 Myongji University, Yongin, Kyonggido 449- 728, Korea 
49 Nagasaki Institute of Applied Science, Nagasaki-shi, Nagasaki 851-0193, Japan 
50 Nara Women’s University, Kita-uoya Nishi-machi Nara 630-8506, Japan 
51 National Research Nuclear University, MEPhI, Moscow Engineering Physics Institute, Moscow, 115409, Russia 
52 University of New Mexico, Albuquerque, New Mexico 87131, USA 
53 New Mexico State University, Las Cruces, New Mexico 88003, USA 
54 Department of Physics and Astronomy, Ohio University, Athens, Ohio 45701, USA 
55 Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA 
56 IPN-Orsay, Univ. Paris-Sud, CNRS/IN2P3, Universite Paris-Saclay, BP1, F-91406, Orsay, France 
57 Peking University, Beijing 100871, People’s Republic of China 
58 PNPI, Petersburg Nuclear Physics Institute, Gatchina, Leningrad region, 188300, Russia 
59 RIKEN Nishina Center for Accelerator-Based Science, Wako, Saitama 351-0198, Japan 
60 RIKEN BNL Research Center, Brookhaven National Laboratory, Upton, New York 11973-5000, USA 
61 Physics Department, Rikkyo University, 3-34-1 Nishi-Ikebukuro, Toshima, Tokyo 171-8501, Japan 
62 Saint Petersburg State Polytechnic University, St. Petersburg, 195251 Russia 
63 Universidade de Sao Paulo, Instituto de Fisica, Caixa Postal 66318, Sao Paulo CEP05315-970, Brazil 
64 Department of Physics and Astronomy, Seoul National University, Seoul 151-742, Korea 
65 Chemistry Department, Stony Brook University, SUNY, Stony Brook, New York 11794-3400, USA 
66 Department of Physics and Astronomy, Stony Brook University, SUNY, Stony Brook, New York 11794-3800, USA 

67 University of Tennessee, Knoxville, Tennessee 37996, USA 
68 Department of Physics, Tokyo Institute of Technology, Oh-okayama, Meguro, Tokyo 152-8551, Japan 
69 Center for Integrated Research in Fundamental Science and Engineering, University of Tsukuba, Tsukuba, Ibaraki 305, Japan 

70 Vanderbilt University, Nashville, Tennessee 37235, USA 
71 Waseda University, Advanced Research Institute for Science and 
Engineering, 17 Kikui-cho, Shinjuku-ku, Tokyo 162-0044 > Japan 
' 2 Weizmann Institute, Rehovot 76100, Israel 

' 3 Institute for Particle and Nuclear Physics, Wigner Research Centre for Physics, Hungarian 
Academy of Sciences (Wigner RCP, RMKI) H-1525 Budapest 114, POBox 49, Budapest, Hungary 
74 Yonsei University, IPAP, Seoul 120-749, Korea 

75 University of Zagreb, Faculty of Science, Department of Physics, Bijenicka 32, HR-10002 Zagreb, Croatia 

(Dated: January 26, 2016) 

We report the measurement of cumulants ( C n , n= 1... 4) of the net-charge distributions measured 
within pseudorapidity (|»j| < 0.35) in Au+Au collisions at ,ys NN =7.7-200 GeV with the PHENIX 
experiment at the Relativistic Heavy Ion Collider. The ratios of cumulants (e.g. C 1 /C 2 , C 3 /C 1 ) 
of the net-charge distributions, which can be related to volume independent susceptibility ratios, 
are studied as a function of centrality and energy. These quantities are important to understand 
the quantum-chromodynamics phase diagram and possible existence of a critical end point. The 
measured values are very well described by expectation from negative binomial distributions. We 
do not observe any nonmonotonic behavior in the ratios of the cumulants as a function of collision 
energy. The measured values of C 1 /C 2 = p/<r 2 and C 3 /C 1 = Sa 3 /p can be directly compared to 
lattice quantum-chromodynamics calculations and thus allow extraction of both the chemical freeze- 
out temperature and the baryon chemical potential at each center-of-mass energy. The extracted 
baryon chemical potentials are in excellent agreement with a thermal-statistical analysis model. 

PACS numbers: 25.75.Dw 


One of the main goals in the study of relativistic heavy 
ion collisions is to map the quantum chromodynamics 


(QCD) phase diagram at finite temperature (T) and 
baryon chemical potential ( pb ) [!]■ Although the ex- 
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act nature of the phase transition at finite baryon den¬ 
sity is still not well established, several models suggest 
that, at large /is and low T, the phase transition be¬ 
tween the hadronic phase and the quark-gluon-plasma 
(QGP) phase is of first order US and that at high T 
and low /is there is a simple cross over from the QGP to 
hadronic phase The point at which the first-order 

phase transition ends in the T — ^lb plane is called the 
QCD critical end point (CEP), which is one of the cen¬ 
tral targets of the Relativistic Heavy Ion Collider (RHIC) 
beam-energy-scan program. Several calculations also re¬ 
ported the possible existence of the CEP in the T — hb 
phase diagram ESS- 

RHIC at Brookhaven National Laboratory has pro¬ 
vided a large amount of data from Au+Au collisions 
at different colliding energies, which gives us a unique 
opportunity to scan the T — fis plane and investigate 
the possible existence and location of the CEP. In the 
thermodynamic limit, the correlation length (£) diverges 
at the CEP [lj. Event-by-event fluctuations of various 
conserved quantities, such as net-baryon number, net- 
charge, and net-strangeness are proposed as possible sig¬ 
natures of the existence of the CEP loj- 13 ■ It has been 
shown in lattice QCD that with a next-to-leading-order 
Taylor series expansion around vanishing chemical po¬ 
tentials, the cumulants of charge-fluctuations are sensi¬ 
tive indicators for the occurrence of a transition from the 
hadronic to QGP phase MM- Typically, the variances 
of net-baryon, net-charge, and net-strangeness distribu¬ 
tions are proportional to £ as o 2 (=C 2 )={(SN) 2 ) ~ £ 2 ;9i], 
where N is the multiplicity, SN = N — fi and =G \) is 
the mean of the distribution. 

Recent calculations reveal that higher cumulants of 
the fluctuations are much more sensitive to the proxim¬ 
ity of the CEP than earlier measurements using second 
cumulants (a 2 ) 0,0. The skewness ( S ) and kurto- 
sis (ft) are related to the third and fourth moments S 
(= C 3 /C 2 3/2 ) = (( 5N) 3 )/o 3 ~ ? 4 - 5 and k(=C 4 /C 2 ) = 
((6N) 4 )/o 4 — 3 ~ £ 7 . The ratio of the various order 
(n) of cumulants (C n ) and conventional values (//, o, S 
and ft) can be related as follows: \x/o 2 = Ci/C 2 , So 
= C 3 /C 2 , no 2 = C 4 /C 2 , and So 3 /n = C 3 /Ci- Be¬ 
cause £ diverges at the CEP, the ratios of cumulants 
So and kct 2 should rise rapidly when approaching the 
CEP H1E3. The cumulants of conserved quantities 
of net-baryon, net-charge, and net-strangeness obtained 
from lattice QCD calculatio ns fl3 . Mm and a hadron 
resonance gas (HRG) model [18( are related to the gener¬ 
alized susceptibilities of n-th order (y n ) associated with 
the conserved quantum numbers as n/ 0 2 ~ X^Vx^j 
So ~ X^ 3 Vx^ 2 \ So 3 /fi ~ X^ 3 Vx^\ an< i K(j2 ~ X (4) /X (2) - 
One advantage of measuring n/ 0 2 , So, So 3 //i, and no 2 
is that the volume dependence of p,, o, S, and ft cancel 
out in the ratios, hence theoretical calculations can be 
directly compared with the experimental measurements. 
These cumulant ratios can also be used to extract the 


freeze-out parameters and the location of the CEP [14| . 
Net-electric charge fluctuations are more straightforward 
to measure experimentally than net-baryon number fluc¬ 
tuations, which are partially accessible via net-proton 
measurement [lj|. While net-charge fluctuations are not 
as sensitive as net-baryon fluctuations to the theoretical 
parameters, both measurements are desirable for a full 
understanding of the theory. 

We report here precise measurements of the energy and 
centrality dependence of higher cumulants of net-charge 
multiplicity (AN c4l = N + — N~) distributions measured 
by the PHENIX experiment at RHIC in Au+Au colli¬ 
sions at y/s NN = 7.7, 19.6, 27, 39, 62.4, and 200 GeV. 
These measurements cover a broad range of /is in the 
QCD phase diagram. 

The PHENIX detector is composed of two central spec¬ 
trometer arms, two forward muon arms, and global de¬ 
tectors [2(ij |. In this analysis, we use the central arm spec¬ 
trometers, which cover a pseudorapidity range of \r]\ < 
0.35. Each of the two arms subtends tt/2 radians in az¬ 
imuth and is designed to detect charged hadrons, elec¬ 
trons, and photons. For data taken at y/s NN = 62.4 
and 200 GeV in 2010 and 2007, respectively, the event 
centrality is determined using total charge deposited in 
the beam-beam counters (BBC), which are also used for 
triggering and vertex determination. For lower energies 
(y/^NN = 39 GeV and below) the acceptance of the BBCs 
(3.0 < I 77 I < 3.9) are within the fragmentation region, so 
alternate detectors must be employed. For data taken 
y/$NN = 39 and 7.7 GeV in 2010, centrality is de¬ 
termined using the total charge deposited in the outer 
ring of the reaction plane detector (RXNP), which cov¬ 
ers 1.0 < \r]\ < 1.5 21|. For data taken at = 19.6 

and 27 GeV in 2011, the RXNP was absent, so centrality 
is determined using the total energy of electromagnetic 
calorimeter (EMCal) clusters to minimize the correlation 
with the charge of the tracks measured in the same ac¬ 
ceptance. More details on the procedure are given in [221. 
The analyzed events for the above mentioned energies are 
within a collision vertex of |Z verte x| < 30 cm. The num¬ 
ber of analyzed events are 2M, 6 M, 21M, 154M, 474M, 
and 1681M for y/s NN = 7.7, 19.6, 27, 39, 62.4, and 200 
GeV Au+Au collisions, respectively. 

The number of positively charged (N + ) and negatively 
charged ( N ~) particles measured on an event-by-event 
basis are used to calculate the net-charge (AV c h) dis¬ 
tributions for each collision centrality and energy. The 
charged-particle trajectories are reconstructed using in¬ 
formation from the drift chamber and pad chambers 
(PCI and PC3). A combination of reconstructed drift- 
chamber tracks and matching hits in PCI are used to 
determine the momentum and charge of the particle. 
Tracks having a transverse momentum ( pr ) between 0.3 
and 2.0 GeV/c are selected for this analysis. The ring 
imaging Cerenkov detector is used to reduce the electron 
background resulting from conversion photons. To fur- 
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ther reduce the background, selected tracks are required 
to lie within a 2.5er matching window between track pro¬ 
jections and PC3 hits, and a 3<r matching window for the 
EMCal. 
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FIG. 1. (Color online). Uncorrected net-charge (A7V c h) distri¬ 
butions, within \r/\ < 0.35 for different energies, from Au+Au 
collisions for (a) central (0%-5%) and (b) peripheral (55%— 
60%) centrality, (c)-(f) are the efficiency corrected cumu- 
lants of net-charge distributions as a function of (N par t) from 
Au+Au collisions at different collision energies. Systematic 
uncertainties on moments are shown for central (0%-5%) col¬ 
lisions. 


Figure Oja) an( i (b) show A7V c h distributions in 
Au+Au collisions for central (0%-5%) and periph¬ 
eral (55%-60%) collisions at different collision energies. 
These AIV c h distributions are not corrected for recon¬ 
struction efficiency. The centrality classes associated 
with the average number of participants ({N palt )) are de¬ 
fined for each 5% centrality bin. These classes are deter¬ 
mined using a Monte-Carlo simulation based on Glauber 
model calculations with the BBC, RXNP, and EMCal 
detector response taken into account 22j, 123J . 

The AN ch distributions are characterized by cumu- 
lants and related quantities such as p, a, S, and k, which 
are calculated from the distributions. The statistical un¬ 
certainties for the cumulants are calculated using the 
bootstrap method 24]. Corrections are then made for 
the reconstruction efficiency, which is estimated for each 
centrality and energy using the hijing1.37 event gener¬ 
ator [2f| and then processed through a GEANT simula¬ 
tion with the PHENIX detector setup. For all collision 


energies, the average efficiency for detecting the parti¬ 
cles within the acceptance varies between 65%-72% and 
76%-85% for central (0%-5%) and peripheral (55%-60%) 
events, respectively with 4%-5% variation as a function 
of energy. The efficiency correction applied to the cu¬ 
mulants is based on a binomial probability distribution 
for the reconstruction efficiency [26J. The efficiency cor¬ 
rected p, a, S, and k as a function of (A r part ) are shown 
in panels (c-f) of Fig. [1] 

The p and a for net-charge distributions increase with 
increasing (A T part ), while S and k decrease with increas¬ 
ing (Ap ar t) for all collision energies. At a given (A part ) 
value, p, S, and k of net-charge distributions decrease 
with increasing collision energy. However, the width (a) 
of net-charge distributions increases with increasing col¬ 
lision energy indicating the increase of fluctuations in the 
system at higher y/s NN . 
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FIG. 2. (Color online). (Apart) dependence of efficiency cor¬ 
rected (a) p/a 2 , (b) Sa, (c) na 2 , and (d) Sir 3 / p of net-charge 
distributions for Au+Au collisions at different collision ener¬ 
gies. Statistical errors are shown along with the data points 
while systematic uncertainties are shown for (0%-5%) colli¬ 
sions. 


The systematic uncertainties are estimated by: (1) 
varying the Z velte x cut to less than ±10 cm; (2) varying 
the matching parameters of PC3 hits and EMCal clus¬ 
ters with the projected tracks to study the effect of back¬ 
ground tracks originating from secondary interactions or 
from ghost tracks; (3) varying the centrality bin width to 
study nondynamical contributions to the net-charge fluc¬ 
tuations due to the finite width of the centrality bins 27- 
Efij ; and (4) varying the lower pr cut. The total system¬ 
atic uncertainties estimated for various cumulants for all 
energies are: 10%-24% for p, 5%-10% for a, 25%-30% 
for S, and 12%-19% for k. The systematic uncertainties 
are similar for all centralities at a given energy and are 
treated as uncorrelated as a function of yjs NN . For clar¬ 
ity of presentation, the systematic uncertainties are only 
shown for central (0%-5%) collisions. 

Figure [2] shows the (A part ) dependence of p/a 2 , Sa, 
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FIG. 3. (Color online). The energy dependence of efficiency 
corrected (a) /i/a 2 , (b) Sa, (c) kct 2 , and (d) Sa 3 / n of net- 
charge distributions for central (0%-5%) Au+Au collisions. 
The error bars are statistical and caps are systematic uncer¬ 
tainties. The triangle symbol shows the corresponding effi¬ 
ciency corrected cumulant ratios for net-charge, from NBD 
fits to the individual N + and N~ distributions. 


na 2 , and Sa 3 / fi{= (S'cr)/(/x/cr 2 )) extracted from the 
net-cliarge distributions in Au+Au collisions at differ¬ 
ent y/s NN . The results are corrected for the reconstruc¬ 
tion efficiencies. Statistical uncertainties are shown along 
with the data points. The systematic uncertainties are 
constant fractional errors for all centralities at a particu¬ 
lar energy, hence they are presented for the central ( 0 %- 
5%) collision data point only. The systematic uncertain¬ 
ties on these ratios across different energies varies as fol¬ 
lows: 20%-30% for /i/tr 2 , 15%-34% for Sa, 12%-22% for 
na 2 , and 17%-32% for Sa 3 / fi. It is observed in Fig. [2] 
that the ratios of the cumulants are weakly dependent on 
(Apart) f° r eac h collision energy; the values of \x/a 2 and 
Sa decrease from lower to higher collision energies, while 
the na 2 and Sa 3 /fi values are constant as a function of 
y/s NN within systematic uncertainties. 

The collision energy dependence of li/a 2 , Sa, na 2 and 
Sa 3 /n of the net-charge distributions for central (0%- 
5%) Au+Au collisions are shown in Fig. [3] The statis¬ 
tical and systematic uncertainties are shown along with 
the data points. The experimental data are compared 
with negative-binomial-distribution (NBD) expectations, 
which are calculated by computing the efficiency cor¬ 
rected cumulants for the measured A + and N~ distri¬ 
butions fit with NBD’s respectively, which also describe 
total charge (A + + N~) distributions very well [27, 281] . 
The various order (n = 1, 2, 3 and 4) of net-charge cu¬ 
mulants from NBD are given as C n (AN c h) = C n (N + ) + 
(— l) n C n (N~), where C n (N + ) and C n (N~) are cumu¬ 
lants of N + and A~ distributions, respectively fffil [31], 

The i-i/a 2 and Sa values in Fig. [3][a) and Fig. [3][b), 


respectively both decrease with increasing y/s NN . The 
NBD expectation agrees well with the data. The na 2 
values in Fig. He) remain constant and positive, between 
1.0 < KtJ 2 < 2.0 at all the collision energies within the 
statistical and systematic uncertainties. However, there 
is ~ 25% increase of na 2 values at lower energies com¬ 
pared to higher energies above y/s NN = 39 GeV, which 
is within the systematic uncertainties. These data are 
in agreement with a previous measurement [32], but pro¬ 
vide a more precise determination of the higher cumulant 
ratios, verified by the NBD method of correcting for ef¬ 
ficiency, which is simple and analytical for all cumulant 
ratios with the standard binomial correction [26|. The 
Sa 3 /\i values in Fig.[3][d) remain constant at all collision 
energies within the uncertainties and are well described 
by the NBD expectation. From the energy dependence 
of [i/a 2 , Sa, na 2 , and Sa 3 /fi, no obvious nonmonotonic 
behavior is observed. Although both previous measure¬ 
ments by STAR [§3, El use the pseudorapidity range 
1 77 1 < 0.5, compared to the present measurement span¬ 
ning |?y| < 0.35, these measurements are all within the 
central rapidity region and are expected to be valid for 
comparison to lattice QCD calculations. The efficiency 
corrected results for the cumulant ratios n/a 2 , Sa, and 
na 2 remain the same within statistics whether each sin¬ 
gle arm of the PHENIX central spectrometer (azimuthal 
aperture 5<j> = 7 r/ 2 ) or both arms {5cj> = n) are used. This 
is a clear verification of the insensitivity of measured cu¬ 
mulant ratios to volume effects. 
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FIG. 4. (Color online). The energy dependence of the 
chemical freeze-out parameter /r_e. The dashed line is the 
parametrization given in Ref. [3J] and the other experimental 
data are from Ref. j34j] and references therein. 


The precise measurement of both [i/a 2 and Sa 3 /fi 
in the present study allow both /is and T/ to be de¬ 
termined, unlike a previous calculation in Ref. 35. .371], 
which was only able to use the [i/a 2 measurement from 
Ref. [32;]. The comparison of Sa 3 /n for different y/s//// 
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TABLE I. Freeze-out Tf and ps vs. y/s NN in the range 27 < y/s NN < 200 GeV from this work compared to ps values from 
Ref. [35], which used STAR net-charge cumulant measurements from Ref. [32;] for ps; with 140 MeV < Tf < 150 MeV obtained 
from STAR net-proton measurement in Ref. [33] by averaging Sa 3 /fi over y/s NN = 27, 39, 62.4 and 200 GeV. 



PHENIX + Ref. [14, 36] 

PHENIX + Ref. [37] 

STAR ± Ref. [35] 

(GeV) 

T S (MeV) 

ps (MeV) 

Tf (MeV) 

ps (MeV) 

ps (MeV) 

27 

164 ± 6 

181 ±21 

160 ±6 

184 ±21 

136 ± 13.8 

39 

158 ±5 

114 ± 13 

156 ±5 

118 ± 10 

101 ± 10 

62.4 

163 ±5 

71 ±8 

159 ±5 

74 ±8 

66.6 ± 7.9 

200 

163 ±8 

27 ±5 

159 ±8 

25 ±7 

22.8 ±2.6 


with the lattice calculations (Fig. 3(b) in Ref. 14, 36]) 
enables us to extract the chemical freeze-out temperature 
(Tf). Furthermore, ps can be extracted by comparing 
the measured p/er 2 ratios with the lattice calculations of 
R 12 = p/(x 2 (Fig.3(a) in Ref. 14, [36}). The extracted 
Tf and ps values are listed in Tabic |T] The Tf and p^ 
extracted using the lattice calculations in the continuum 
limit from Ref. [37| are also depicted in Table [I] The ex¬ 
tracted freeze-out parameters using different lattice re¬ 
sults agree very well. However, the extracted Tf are 
2-4 MeV lower using Ref. 371 than with Ref. 14, [3(1, 
which is well within the stated uncertainties. The de¬ 
tailed freeze-out parameter extraction procedure is given 


in Ref. [1J, 35, 37]. This is a direct combination of ex¬ 


perimental data and lattice calculations to extract phys¬ 
ical quantities. The y/s NN dependence of ps shown in 
Fig. [4] is in agreement with the thermal-statistical anal¬ 
ysis model of identified particle yields [34| . The ps ex¬ 
tracted in the present net-charge measurement and the 
values reported in [35] are in agreement within stated un¬ 
certainties, with some tension at y/s NN = 27 GeV. Avail¬ 
able lattice results allow extraction of ps and Tf from 
= 27 GeV and higher using the present net-charge 
experimental data. Other recent calculations [H, [ 3 ^ ] 
have used both net-proton and net-charge measurements 
to estimate the freeze-out parameters. 


In summary, fluctuations of net-charge distributions 
have been studied using higher cumulants (p, o\ S, and 
At) for \r]\ < 0.35 with the PHENIX experiment in Au+Au 
collisions ranging from y/s^ = 7.7 to 200 GeV. The ra¬ 
tios of cumulants (p/cr 2 , Sa, Atcr 2 , and Sa 3 /fT) have been 
derived from the individual cumulants of the distribu¬ 
tions studied as a function of (7V part ) and y/s NN . The 
/la,/ct 2 and Sa values decrease with increasing collision en¬ 
ergy and are weakly dependent on centrality, whereas Atcr 2 
and Sa 3 /fj, values remain constant over all collision ener¬ 
gies within uncertainties. The efficiency corrected values 
from the NBD expectation reproduce the experimental 
data. These data are in agreement with a previous mea¬ 
surement [ 32 }, but provide more precise determination of 
the higher cumulant ratios Sa and a ta 2 . In the present 
study we do not observe any significant nonmonotonic 
behavior of p/cr 2 , Sa, na 2 , and So 3 j[i as a function of 


collision energies. Comparison of the present measure¬ 
ments together with the lattice calculations enables us to 
extract the freeze-out temperature Tf and baryon chem¬ 
ical potential (hb) over a range of collision energies. The 
extracted /is values are in excellent agreement with the 
thermal-statistical analysis model 341. 
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